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ABSTRACT 

Some hot, massive, population-I Wolf-Rayet (WR) stars of the carbon subclass are known to 
be prolific dust-producers. How dust can form in such a hostile environment remains a mystery. 
Here we report the discovery of a relatively cool, extended, multi-arc dust envelope around the 
star WR112, most likely formed by wind-wind collision in a long-period binary system. We 
derive the binary orbital parameters, the dust temperature and the dust mass distributions in 
the envelope. We find that amorphous carbon is a main constituent of the dust, in agreement 
with earlier estimates and theoretical predictions. However, the characteristic size of the dust 
grains is estimated to be ~ 1 /loti, significantly larger than theoretical limits. The dust production 
rate is 6.1 x lO~'M 0 yr _1 and the total detectable dust mass is found to be about 2.8 x 1O” 5 M 0 
(for d = 4.15 kpc). We also show that, despite the hostile environment, at least ~20% of the 
initially-formed dust may reach the interstellar medium. 

Subject headings: infrared: stars — stars: Wolf-Rayet — stars: winds, outflows 


1. Introduction 

Classes of objects known to produce significant 
quantities of dust include asymptotic giant-branch 
stars, red giants, novae, supernovae (each of which 
contributes about equal rates, ~ lO _3 M 0 yr _1 , 
of dust in the whole Galaxy), along with plane¬ 
tary nebulae and proto-stars (each of which yields 
~10 times less; Dwek 1985). Surprisingly, some 
hot, massive, population I WR stars also form 
dust (van der Hucht, Williams, & The et al. 
1987). WR stars are the evolved descendants 
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of massive O-type stars. They consist mainly 
of He-burning cores surrounded by hot envelopes 
that drive fast, dense winds with average mass- 
loss rates ~ 10 -5 M 0 yr _1 and terminal velocities 
Vqo ~ 1000 — 4000 kms^ 1 . There are three suc¬ 
cessive WR phases, WN, WC and WO, charac¬ 
terized by the dominant emission lines of N, C, 
and O, respectively, in their optical spectra. All 
dust-making WR stars belong to the carbon-rich, 
hydrogen-poor WC subclass (Williams 1995). 
The WR dust-makers are remarkable for two main 
reasons: (1) the absolute rate of formation is very 
high, up to 10 -6 M 0 yr _1 in dust alone, and occur¬ 
ring in either a periodic (due to enhanced wind- 
wind compression at or near periastron passage in 
long-period WC + O binaries with eccentirc or¬ 
bits) or sustained fashion (in single late-type WC 
stars or moderately short-period WC + O binaries 
with circular orbits); and (2) the dust is formed in 
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a hot, extremely hostile environment, where the 
formation process is still unknown. 

Here we report on near and mid-infrared imag¬ 
ing observations of the dust envelope surrounding 
the star WR 112 (spectral class WC9). The mor¬ 
phology of this envelope provides clues to the na¬ 
ture of the stellar system, while the photometry 
allows us to estimate the properties of the dust 
and the total dust mass. 

2. Observations 

WR112 belongs to the group of 5 WC stars with 
the densest known dust envelopes (van der Hucht 
et al. 1996). We observed WR112 along with 
WR104 and WR118 with the University of Florida 
mid-IR imager ’’OSCIR” at the Gemini-North 8m 
telescope using the medium-band 7.9/xm, 12.5/xm 
and 18.2 fxm filters on 7 May, 2001. For WR112 
we supplemented the mid-IR data with near-IR 
narrow-band h and k images taken at the CFHT 
on 27 June 1999 with the adaptive optics bonnette 
and KIR camera, and broad-band L ' images ob¬ 
tained at the IRTF on 28 May 2000 in the ‘movie 
burst’ (2 x 8192 0.01-sec exposures) mode of the 
NSFCAM imager. We also acquired images of 
point sources (PSF reference stars), usually im¬ 
mediately before and after the principal target, 
for image restoration. The original data cubes 
were split into individual frames, quickly checked 
and selected for image quality using basic image 
statistics for the OSCIR and CFHT data (initially 
compensated for atmospheric turbulence), or care¬ 
fully selected using a specifically designed statis¬ 
tical procedure for the NSFCAM images (taken 
without compensation). The individual images 
were then appropriately shifted to form a final im¬ 
age for each wavelengtlr/target/instrument. The 
final images were then restored using the images 
of the PSF reference stars and a maximum en¬ 
tropy algorithm (the ‘stsdas.analysis.restore.mem’ 
task of IRAF 4 for the CFHT data, enabling us 
to reach 1 pixel (0.018") resolution, or the ‘sts¬ 
das.analysis.restore.lucy’ task of IRAF for the OS¬ 
CIR and NSFCAM images, thus reaching 1-pixel 
resolution of 0.089" and 0.055", respectively). 

4 IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agree¬ 
ment with the National Science Foundation. 


3. Dust properties 

While WR104 and WR118 show only modest 
extensions on scales of several 0.1", combination 
of the original mid-IR images of WR112 reveals 
three sets of regularly spaced arc-like structures 
extending out 3" from the central star (Fig. la). 
It is apparent that the dust temperature gradu¬ 
ally falls off outwards. The central source and 
two closest arcs form a W-like configuration, the 
middle part of which is also seen in the near-IR 
(Fig. lb). The near-IR size of the central source, 
FWHM=0.089"±0.004", is consistent with earlier 
measurements (Ragland & Richichi 1998; Mon- 
nier et al. 2002): FWHM=0.059" - 0.073". The 
near-IR SW elongation probably corresponds to 
the previously detected (Ragland & Richichi 1998) 
small-scale asymmetry in the dust envelope. 

WR112 is the fourth case in which dust around 
a WR star has been mapped in some detail, with 
all the previous examples occurring in WC+O bi¬ 
naries (Marchenko, Moffat, & Grosdidier 1999; 
Monnier, Tuthill, & Danchi 1999; Tuthill, Mon- 
nier & Danchi 1999). Note that 6 other WR stars 
have been shown to posess an extended (presum¬ 
ably spherically-symmetric) envelope (Yudin et al. 
2001; Monnier et al. 2002). WR112 is also known 
to be a variable non-thermal radio emission source 
(Chapman et al. 1999; Monnier et al. 2002), which 
might indicate binarity (Dougherty & Williams 
2000). There are few viable scenarios for dust pro¬ 
duction in WR112. For example, the primordial 
circumstellar envelope would not have survived in 
the hard UV radiation-field of the WR progeni¬ 
tor. Periodic outbursts driven by stellar oscilla¬ 
tions cannot explain the observed dust envelope- 
shape without invoking a factor of two azimuthal 
variation in the terminal-wind velocity along with 
an extreme latitudinal dependence of the mass- 
loss. Therefore, we assume that the dust is formed 
in the wind-wind collision zone of a massive binary 
system (thus resembling a collimated but rotating 
‘beam’ for a distant observer) and expelled radi¬ 
ally outward at a constant, terminal-wind veloc¬ 
ity i>oo = 1200 fcms -1 (Rochowicz & Niedzielski 
1995). We then calculate the resulting locus of 
maximum dust density and fit it to the observed 
dust distribution in the restored mid-IR images, 
taking the necessary orbital elements, P, e, i, u>, f 1 
as free parameters. We find a family of minimum- 
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X 2 solutions for counter-clockwise orbits, in which 
coupling of e-i and w-Sl leads to some degree of 
degeneracy. Instead of relying solely on the sin¬ 
gle formal best fit, we therefore construct a syn¬ 
thetic solution by averaging all the possible com¬ 
binations of solutions for which the % 2 values lie 
within 5% of the absolute minimum. This average 
solution has period P = 24.8 ± 1.5 yr, eccentricity 
e = O.lliO.ll, orbital inclination i = 38.0 ±3.8 , 

O O 

periastron angle to = 269.6 ±18.5 and orbital ori¬ 
entation O = 33.5 ± 7.9 (1 -a errors). In Fig. 2 
we plot the solution with the highest possible, but 
still plausible, eccentricity: e = 0.40. This par¬ 
ticular model yields a dust locus that is actually 
very close to that derived from the synthetic so¬ 
lution. In general, the model provides fairly good 
fits to the SW arms, although it fails to explain the 
pronounced NE excursions, as well as the bright 
SW protuberance. The latter extends in the di¬ 
rection toward a stellar visual companion ~ 1" 
from WR112 (Wallace, Moffat & Shara 2002, and 
our Fig.2), which is not seen at near-IR, and is 
therefore probably also invisible at mid-IR wave¬ 
lengths. Clearly the dust-production rate depends 
on orbital phase, diminishing as expected around 
apastron but surprisingly also absent around pe¬ 
riastron, contrary to the situation in the eccen¬ 
tric, long-period colliding-wind WR±0 binaries 
WR137 and WR140 (van der Hucht, Williams, 
& Morris 2001; Marchenko, Moffat, & Grosdidier 
1999). The average thickness, Sr, of the dust arcs 
in the deconvolved images relative to the separa¬ 
tion between the successive arcs, Ar (measured 
at the same orbital phase), gives an estimate of 
the shock-cone opening angle: 0 ~ 27 t x Sr/Ar ~ 

O O 

60 —110 , similar to that found in other colliding- 
wind WR±0 systems (Marchenko et al. 1997; 
Bartzakos, Moffat, & Niemela 2001). In Table 1 
we list the calculated orbital elements along with 
the main characteristics of the dust envelope and 
the assumed parameters. 

To derive the basic dust properties, which for¬ 
tunately do not depend sensitively on the adopted 
orbital parameters, we first deproject all apparent 
distances by applying the synthetic average orbit- 
solution. We calibrate the restored images us¬ 
ing all available IR-flux measurements of WR112 
(Williams, van der Hucht, & The 1987; van der 
Hucht et al. 1996), assuming that the system has 
not experienced any large-amplitude IR outbursts. 


In the mid-IR range, the relatively large apertures 
used in previous ground- and space-based obser¬ 
vations will tend to smooth out the variability by 
averaging over the output of many dust-formation 
cycles. We deredden the near-IR data by apply¬ 
ing a Galactic interstellar-extinction model (Are- 
nou, Grenon, & Gomez 1992), then calculate the 
net dust fluxes after subtracting the stellar compo¬ 
nent, approximated as F\ ~ A _a , with a = 2.97 
(Morris et al. 1993). Then we use the relation 
(Hildebrand 1983): 


M d (r) 


ApF{\, r)d 2 a 
3B(\, T(r))Q e (X, a) ’ 


(1) 


where M d is the dust mass, F( A, r) is the measured 
flux at a given position r, p is the grain density, 
d=4.15 kpc is the distance to WR112 (van der 
Hucht 2001), a is the grain size, B(X,T(r)) is the 
black-body emissivity, and Q e (X,a) is the grain 
emission coefficient. We adopt the general relation 
Q e (X,a)/a = kX~ s along with a T = Tor -/3 tem¬ 
perature profile (Williams, van der Hucht, & The 
1987) and measure the fluxes F(X,r) over regions 
of fixed 7x7 pixel size (0.623" x 0.623"), assigning 
15% /30% errors to F(X,r) in the inner/outer¬ 
most arcs. Taking the ratios of fluxes at every 
given position in the envelope (except the very 
central part, where the spectral distribution can¬ 
not be represented by a single temperature) in or¬ 
der to reduce the number of free parameters and 
to be impervious to large positional fluctuations 
in F(X,r), we obtain via ^-minimization the fol¬ 
lowing estimates: s = 0.90+o'Ji, T 0 = 320^ K, 
/3 = 0.40io!o5 (To at 1"; 95% confidence intervals). 
The dust emissivity index s, though loosely con¬ 
strained, is close to the range expected for amor¬ 
phous carbon in the near-mid-IR domain (Borgh- 
esi, Bussoletti, & Colangeli 1985; Suh 2000). On 
the other hand, the absence of the telltale 9.7 pm 
silicate emission feature in the WR112 spectrum 
(van der Hucht et al. 1996) rules out the pres¬ 
ence of any silicate-based mixture. The derived 
temperature profile follows the (3 = 0.4 depen¬ 
dence expected in the case of thermal equilibrium. 
Adopting p = 2.0 g cm~ 3 for the amorphous car¬ 
bon dust, we estimate the dust mass separately 
for each 7x7 pixel 2 area of each of the 3 mid-IR 
images, then average the results (Fig. 3). We find 
a total M d = 5.5 x 10 28 g, integrated over the ob¬ 
served envelope. To compare our result with an 
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earlier theoretical estimate of the mass we (a) ac¬ 
count for the difference in the adopted distances 
to WR112 and (b) assume that the observed dust 
is created for 15 yrs (see Fig. 2) during each 25-yr 
orbital cycle, with 3 complete cycles accounted for. 
A model (Zubko 1998) with the assumption of a 
spherically-symmetric envelope gives M ( j = 3.8 x 
lO~'M 0 yr _1 (converted to d=4.15 kpc), while 
our estimate is M,i = 6.1 x 10“' M 0 yr -1 . Calcu¬ 
lation of the integrated mass-ratios, M(innermost 
region+arcs) : M(middle arcs) : M(outer-most 
arcs)=1.0 : 0.6 : 0.2, shows that the dust is gradu¬ 
ally destroyed while leaving the system. Both the 
intense UV radiation field and sputtering caused 
by the high drift-velocities of the dust grains 
(Zubko 1998) may be responsible for the dust de¬ 
struction. Nevertheless, the flattening of the trend 
at large radii in Fig. 3 suggests that some dust 
does reach the ISM. 

Finally, we use eq. 1 to evaluate the characteris¬ 
tic size of the dust grains, by comparing via least- 
square minimization, the observed fluxes to the 
prediction of a simple model assuming that the en¬ 
velope is populated by similar-sized spherical par¬ 
ticles. We calculate the Mie absorption coefficients 
(Bohren & Huffman 1998), Q a bs(X,a) = Q e (X,a), 
using the optical constants for amorphous carbon 
(Zubko et al. 1996). As a second free parameter we 
take the integrated (innermost region+arcs) dust 
mass, allowing it to follow the general trend seen 
in Fig. 3. This gives an estimate of Md within 25% 
of the previously calculated value, thus providing 
a good consistency check. We find a character¬ 
istic grain-size a = O.491q 0 (^?7j (95% confidence 
interval), in line with the estimates obtained for 
WR112 from ISO spectroscopy, a ~ 1 fim (Chair 
& Tielens 2001). The large grain size poses a seri¬ 
ous problem for the theory of grain growth, which 
calls for a ~ O.Ol^m (Zubko 1998). 

4. Conclusions 

The recently acquired high-resolution, high 
signal-to-noise near-mid-IR images of the WR star 
WR112 have enabled us to resolve and, for the first 
time, study in detail the extended dust envelope 
thought to be formed in the wind-wind collision 
zone of the long-period, P « 25 yr, binary. A sim¬ 
ple approach involving a minimum of assumptions 
has allowed us to derive the basic dust properties: 


We find that the dust emissivity approximately 
follows the trend expected for amorphous car¬ 
bon. The radial temperature-distribution shows 
the dust envelope to be in thermal equilibrium. 
The dust production-rate corresponds to ~ 6% of 
the total mass-loss rate, if for the latter we assume 
a typical value of 10 -5 M 0 yr -1 . The characteris¬ 
tic size of the dust particles, a ~ 1/rm, turns out to 
be much larger than expected from state-of-the- 
art models of dust growth. We also show that, 
despite the harsh conditions, ~ 20% of the newly 
formed dust can survive and escape the system, 
thus enriching the interstellar medium. This may 
have a direct implication for the process of heavy- 
element enrichment in the early Universe, when 
massive stars completely dominated the stellar 
scene. 
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Fig. 1.— a: False-color image of WR112 pro¬ 
duced by log-scaling, normalizing and combining 
the 7.9, 12.5 and 18.2/xm images as blue, green 
and red. The faintest details are ~ 10~ 3 of the 
central source. The inset shows a point-spread- 
function standard at the same scale, used for im¬ 
age restoration. Fig. lb: Zoom on the central 
source: combination of the maximum-entropy de¬ 
convolved, log-scaled and normalized images of 
WR112, with the narrow-band images h (Ao = 
1.61/iTO, FWHM=0.02/jm) and k (Ao = 2.17 /im, 
FWHM=0.02/zm) as blue and green, respectively, 
and the broadband L' image (Ao = 3.75/zm) as 
red. The deconvolved and appropriately enlarged 
7.9 fim image is shown in contours. Small black 
stars mark the same central position on each im¬ 
age. 
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Fig. 2.— Deconvolved 18.2m image of WR112 
with the overplotted model fit corresponding to 
the high-eccentricity case, P = 23.5yr, e = 0.40, 

o o o 

i = 35 , = 40 , u> = 255 (for the average de¬ 

rived parameters see Table 1), with the red line 
tracing the locus of visible dust, then gray line 
beyond. The grossly enlarged, projected orbit is 
shown in the lower left corner. The red dotted 
line points toward periastron. A small blue cross 
depicts the position of the visual companion. 


Fig. 3.— The integrated (within 7x7 pixel ar¬ 
eas) dust mass in grams. Filled triangles mark 
the central part of the image along with the two 
innermost arcs. Filled (open) squares/circles cor¬ 
respond to the SW (NE) middle/outer-most arc. 
The dotted line shows the trend expected in the 
case of mass conservation. 
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Table 1 


Orbital parameters (‘synthetic solution’), basic dust characteristics and assumed values. 


Parameter 

Value 

Errors 3- 

p, yr 

24.8 

±1.5 

e 

0.11 

±0.11 


o 

o 

i 

38.0 

±3.8 


o 

O 

UJ 

269.6 

±18.5 


o 

O 

n 

33.5 

±7.9 

S 

0.90 

-0.21, ±0.22 

T 0 ?K 

320 

-13,±12 

0 

0.40 

-0.05, ±0.05 

M d , g 

5.5 X 10 28 

— 

M d , M 0 yr -1 

6.2 x 10~ 7 

— 

a , fim 

0.49 

-0.11, ±0.11 

Vdust! kms -1 

1200 

assumed 

Pdust, gem -3 

2.0 

assumed 

d, kpc 

4.15 

assumed 


Either 1-cr or 95% confidence interval 
















Distance, ” 



